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Abstract 

The deuteron deep inelastic unpolarized structure function F2' is calculated us- 
ing the Wilson operator product expansion method. The long distance behaviour, 
related to the deuteron bound state properties, is evaluated using the Bethe-Salpeter 
equation with one particle on-mass~shell. The calculation of the ratio F'^/2F'^ is 
compared with other convolution models showing important deviations in the re- 
gion of large x. The implications in the evaluation of the neutron structure function 
from combined data on deuterons and protons are discussed. 
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1 Introduction 



Since the discovery of the EMC effect it became clear that the nucleon structure 
functions inside nuclei are different from the ones determined for free nucleons. In the 
quark-parton model point of view, this means that the quark and the anti-quark sea 
distributions, in the nucleon, are changed by the surrounding nucleons. Traditional meth- 
ods @, 1^ that considered only Fermi motion of the nucleons are not sufficient to describe 
the experimental data. Several efforts have been done to explain this deviation, mainly 
along two different approaches. The first is QCD motivated and assumes a rescal- 
ing, i.e., F2{x, = ^2 ^Q^); where ^ is a A-dependent free parameter determined 

by the best fit to the data. The other typical approach 0-[0 based on the nuclear 
degrees of freedom namely, the nucleon-nucleon interaction involving mesonic exchange 
contributions and other nuclear structure effects. In this work special attention is given to 



the method described in Ref . |T3, ^ , which relies on the Wilson operator product expan- 
sion (OPE). This is a systematic and self-consistent approach inspired in nuclear physics 
concepts which will be used here in the calculation of the deuteron structure function F2 ■ 
A comparison [|I^ between the rescaling model and the OPE model through the 
moments M^[F2) of the structure functions, shows that both methods agree in the 
intermediate region of x (0.2 < x < 0.7) corresponding to small values of n. For large n, 
it seems that the OPE is more appropriate since it describes the experimental data in the 
region of large x. The phenomenological nature of the rescaling method results from 
the fact that it is hard to predict the value of within QCD. On the other hand, the OPE 
method has no free parameters and, in particular, the deuteron structure function can be 
obtained in a consistent way. 



Ref. [|I2| provides a general covariant method to describe deep-inelastic scattering 
(DIS) that does not rely on the Bjorken limit. In this model the deuteron structure 
function depends upon the off-shell nucleon structure function and the Bethe-Salpeter 
deuteron wave function. However, there is no factorization, which implies that, in general, 
the convolution formula is not valid. The usual ambiguities related to the off-shellness 
problems in DIS are not present, but the nucleon distribution function, which is an im- 
portant tool for the interpretation of other processes involving the deuteron, is not con- 
sistently obtained within this formalism. 

Some important considerations about parton distribution functions and QCD analysis 
of DIS rely on the experimental structure functions of the nucleons. In general, the 
neutron structure function is extracted from combined experimental data on the proton 



and deuteron |TJ], taking into account only the Fermi motion of the nucleons in the 
deuteron. The question of whether the EMC effect is a feature only of the heavy nuclei 
or should also be considered as leading to an important correction in light nuclei must be 
discussed. In principle one should take into account the binding effects and the meson 
exchange currents, required to describe the EMC effect in heavy nuclei, also in light nuclei. 
It has been shown that the additional corrections are important in the deuteron and 
that the extracted neutron structure function is model dependent due to the nuclear 
contributions. This is particularly important at large x where the discrepancy with the 
data is of the order ~ 25% and exceeds the experimental error bars |T^. It is thus clear 
that, from a nuclear physics point of view, one should consider corrections to the Fermi 
motion of the nucleons even for light nuclei. 

In this paper the deuteron structure function is calculated using a relativistic formal- 
ism, based on the Wilson operator product expansion method, following Ref. M 
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The deuteron Bethe-Salpeter amplitudes are obtained from a quasi-potential equation 
with one particle on-shell . Section ^ is an overview of the formalism and presents 
the relevant expressions used in the calculation. The main features of the relativistic 
deuteron wave function used in our calculations are described in section |^ and in section 
^ we discuss the results of the ratio jlF^ through a comparison with other approaches. 



2 The Formalism 

The basic feature of the Wilson operator product expansion is the possibility of factorizing 
the DIS amplitude into two pieces in the limit ^ cxd, one associated with the long 
distance behaviour and the other with the short distance behaviour. The pertubative 
QCD formalism determines the large behaviour of the amplitude. On the other hand, 
a nuclear physics approach can be used to deal with the target nuclei, where the long 
distance behaviour can in general be predicted from the bound state properties. 

According to the OPE method |T0|, |11| the deuteron structure function F2 satisfies 
the relation 

\" dxx-^F?M'')=Y.c'^aUm^''r . (1) 
•''0 

where x = Q'^/2Po-q is the usual Bjorken scaling variable. The sum in the second member 
of Eq.(|l|) runs over different fields of the theory, according to the twist-2 approximation. 
The Wilson coefficients Cj^^l are target independent, giving the short distance behaviour 
of F^. In the impulse approximation, it is found that these coefficients are identical to 
the moments of the structure function F2 of nucleons (a = A^) or mesons (a = B), 

cSUQ") = Mn{F^) 4'i(Q') = ) . (2) 

Noting that the moments are defined as 

M„(F) = dxx''-^F{x,Q''), (3) 
Jo 

we obtain from Eqs.(|l]),(@) 

M„(Ff ) = M^Fi^) + M„(Ff ) /.f , (4) 

where /i^^ is interpreted as the moment of an effective distribution function of nucleons 
or mesons in the deuteron. In the leading twist approximation, these functions are give 
by 

< Pd I O^^-^" I Pd >= . . P^" (5) 

where O'^^-^" is the set of local operators that provide the basis for the operator prod- 
uct expansion. The matrix elements from Eq.(|^) can be explicitly evaluated using the 



Mandelstam method fT^ and in particular, the following result holds in the impulse 
approximation 



X 



5(pl-ml) ip,.er~' + (1^2) (6) 
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where pi, p2 are the nucleon momenta in the deuteron and k = [pi —p2)/2. The vector 
e is chosen in such a way that = and e (x q where q is the photon momentum with 
q"^ = —Q^. In particular, if the transfered 3-momentum q is along the z-axis one can 
choose e = (1,0,0,-1). This prescription is used in extracting fx^^^ from the twist-2 
operator expansion. 



Expression @ is similar to the one presented in Ref . [|TT| , but includes a 5-function in 
the integrand which means that the distribution function moment fi^^^ is calculated with 
one particle on-mass-shell. The difference is related to the fact that the Bethe-Salpeter 
amplitudes x, x used in this work, are determined by a quasi-potential equation with one 
particle on-mass-shell [1^. This approximation includes some off-shell effects (the other 
particle is completely off-shell) and it reproduces successfully the bound state properties 
including the deuteron form factors ||18| . 



Applying the inverse Mellin transform to Eq. (|J) , defined as 



F{x) = -^ dnx-^M^F), (7) 



and taking into account the explicit expression for n^^^ one recovers the convolution 
formula 

(x) = /' dy f'^{y) F^{x/y) + MEC (8) 



where 



^/^(y) = ^/(0JX(P1,P2)5(PO/5(P1)X(P1,P2)5-^(P2) 

X 5{pI- ml) l^(pi.e) 5 [y - ^) e{-p,.e) 5 [y + 



pi.e 



+ (l<-2) (9) 

is the distribution function of nucleons in the deuteron. The meson exchange currents 
(MEC) also contributes to the deuteron structure function and the corresponding distri- 



bution function can be found in an analogous way |TT 



Notice that f^l^ given by Eq. (|^) takes into account the binding and off-shell effects in 
the deuteron, including the features needed to describe the EMC effect. It can be shown, 
using the Bethe-Salpeter normalization condition, that the distribution function satisfies 
the baryonic number relation, 

dyf^'^{y) = 2. (10) 



The average value of the momentum carried by the nucleons is, 

<y>= C dyyf''/''{y), (11) 
Jo 

which can be written in the form 

<y>=l-5N, (12) 

where 5n is that part of the deuteron momentum carried by the mesons. The value of 
5m controls the magnitude of the binding effects. Since the function f^^^{y) is strongly 
peaked at ?/ = 1/2, in a first approximation one has f^'^{y) ~ 25{y - 1/2), resulting 
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in (5jv = 0. In this case the deuteron is assumed to be made of two non-interacting 
quasi-free nucleons, where binding effects are neglected. This approximation leads to 
F^{x) ~ 2F,^{2x), and corresponds to a non-smearing motion of the nucleons. A better 
estimative of can be obtained if we expand around y =< y > /2in the integrand 
of Eq. (|) . The result is given by 

(x) = 2Fi^{2x/ <y>) + l{<y'>-<y>' /2)^^^^|^ + . . . (13) 



3 Deuteron vertex 

In relativistic field theory we can describe the deuteron by the Bethe-Salpeter (BS) vertex 
Xs , where s = —1, 0, 1 is the deuteron helicity. This function is the fourier transform of 
matrix elements involving the nucleon fermionic fields, 

Xs = J d^xe'''-'' < I T^{x/2) ^(-x/2) | Pd, s > (14) 



and satisfies the homogeneous BS equation |19|. It depends on the deuteron total momen- 



tum Pd and on the nucleons relative momentum k. According to the Gross on-mass-shell 



prescription |T6[ this vertex can be parameterized by, 



Xs = Fi+—{k.O + 

niN 

+ iP2lI^i^iH(+^(k.O), (15) 

niN rriN 

where ^^{Pd, s) are the deuteron polarization functions, a set of spin dependent space-like 
orthogonal vectors satisfying the following identities 



UPd, s)Pl = C,{Pd, s') ^'^{Pn, s) = -S,,^ 

(16) 

Er (Pi., s) C{Pd, s) = -{g^'^ _PDPh, 



m 



D 



The scalar amplitudes F, G, H and Y depend only on the variable a = [{P£,.k)/'m£)f — k"^ 
and are determined covariantly by solving the BS equation with one particle on-mass- 



shell. The solutions can be found in Ref. 20 



The NN interaction used involves the meson exchange (cr, 00,11, p) with the parameter- 



ization taken from the first article in Ref. 1 16]. We also consider two different couplings for 



the pion namely, a pseudovector (PV) and a mixed pseudoscalar-pseudovector (PS-PV) 
coupling, which is characterized by the parameter A G [0, 1] defined in the pion-nucleon 
vertex as 

where is the pion coupling constant. The PV coupling leads to a deuteron wave function 
consisting only of S and D states. The percentage of P states is negligible for small A 
but increases with A and for the PS-PV coupling corresponding to A = 0.4 is 0.45 %. In 
this case the interference of the P and S states can lead to a significant effect. In our 
calculations we consider two solutions, one obtained with A = (BSO) and the other with 
A = 0.4 (BS04). Both solutions are consistent with the deuteron static properties. 
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4 Results and Conclusions 



The extraction of the neutron structure function from combined experimental data on 
the proton and deuteron can now be discussed. The formahsm presented here indicates 
that one should use Eq.(|^)-(P) to analyse this problem consistently in the large region 
where the Wilson operator product expansion is reliable. The on-shell condition given in 
Eq.(ll) corresponds to a well established and successfull procedure in nuclear physics. In 
order to get some physical insight into the corrections introduced by Eq.(y) one has to 
compare the results obtained by different approaches. In Fig. (||) we present the results of 
the ratio /2F^ function of xj^ = (m/)/mjv)x, for the Bethe-Salpeter formalism 
with a simpler deuteron model (BSsm) |]lT|, a calculation using a non-relativistic deuteron 
wave function (NR) fl^, a light-cone calculation (LC) and finally the Bethe-Salpeter 
(BS) formalism with one particle on-mass-shell. In the latter calculations we used the 
convolution formula (|) with the corresponding distribution function given by Eq.(P), 
while in the other approaches the results were obtained using the expansion (13). The 
nucleon structure function F.^ is taken from Ref. ||13| . 

It is clear that for small values of a: at all different approaches agree. However, for 
xn > 0.3 the deviation between the different models increases. The main reason for 
these differences can be attributed to binding effects. In fact, the value of 6iy differs 
substantively in all approaches. It goes from = in LC calculations to 6iy ~ 5.0 x 10^^ 

N 

^ ~ 7.2 X 10-3 and 6n ^ 6 
respectively. 

The deviation between the two curves BSO and BS04 is an effect of the deuteron P 
states, which are present only on the last one. In particular, the BS04 result deviates 
from all the other calculations even for small values of x. This is not surprising since 
none of those models consider the presence of deuteron P states. The present calculations 
indicate that there is an ambiguity in the extraction of the neutron structure function 
from the combined deuteron and proton data. The procedure is model dependent and 
particularly sensitive to the features of the deuteron wave function. 



in NR calculations with the intermediate value S^- ~ 3.9 x 10 ^ for the BSsm. In our 
calculations we obtain 5m ~ 7.2 x 10"^ and 5m — 6.3 x 10"^ for the BS04 and BSO models, 
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Figure 1: Theoretical predictions for the ratio F.^/2F,^ using four different approaches. 
The curves tagged by BSO and BS04 are the results of the present calculation. In the 
first one, we consider a pseudo- vectorial coupling for the pion-nucleon interaction and in 
the second one a mixed coupling (pseudo-vectorial plus pseudo-scalar) as described in 
the text. The remaining curves result from calculations using the simplified deuteron BS 
wave function of Ref.[ll] (BSsm), the non-relativistic deuteron wave function of Ref.[10] 
and the hght-cone approach of Ref.[21]. 
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